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We investigate Rayleigh scattering in dissipative optical lattices. In particular, following recent 
proposals (S. Guibal et al, Phys. Rev. Lett. 78, 4709 (1997); C. Jurczak et al, Phys. Rev. Lett. 77, 
1727 (1996)), we study whether the Rayleigh resonance originates from the diffraction on a density 
grating, and is therefore a probe of transport of atoms in optical lattices. It turns out that this 
is not the case: the Rayleigh line is instead a measure of the cooling rate, while spatial diffusion 
contributes to the scattering spectrum with a much broader resonance. 
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PACS numbers: 42.65.Es,32.80.Pj 



Light scattering y|, i.e. the scattering of photons re- 
sulting from the interaction with a material medium, is 
a technique widely used to determine the properties of 
many different types of media. From the position and 
the width of the scattering resonances it is in fact pos- 
sible to identify the dynamical modes of the system and 
derive the rates of relaxation toward equilibrium. This 
is well exemplified by the Landau-Placzek relation, valid 
for light scattering originating from the density fluctua- 
tions of a medium at thermal equilibrium, which connects 
the strength of the different components of the scattering 
spectrum to the specific heats of the medium at constant 
volume and constant pressure. 

Recently light scattering has been extensively used to 
study the properties of cold atomic samples, and in par- 
ticular it turned out to be an essential tool for the un- 
derstanding of the basic properties of dissipative optical 
lattices Q • The same technique may also apply to far-off- 
resonance non-dissipative optical lattices which are cur- 
rently investigated by many groups in connection with 
Bose-Einstein condensation experiments |3j]. However to 
derive the damping rates of the system from light scatter- 
ing measurements is in general a highly non-trivial task. 
This is especially true for quasi-elastic (Rayleigh) scat- 
tering 0, la, @ j which gives access to the relaxation rates 
of non-propagating material observables. In this work we 
investigate the mechanism behind the Rayleigh scatter- 
ing in dissipative optical lattices, and identify the relax- 
ation process which determines the width of the Rayleigh 
resonance in the scattering spectrum. 

The starting point of the present study is the previ- 
ous claim that Rayleigh resonances may originate from 
the excitation of the atomic density, and consequently 
the width of the Rayleigh line would provide a mea- 
sure of the diffusion coefficients of the atoms in an op- 
tical lattice 0. Following a similar approach, Jurczak 
et al derived values for the diffusion coefficients from 
polarization-selective intensity correlations. 

In our analysis we first assume, on the lines of these 
previous works, that the material observable excited in 
the pump-probe spectroscopy is the atomic density, and 
derive the expected relation between the width of the 
Rayleigh resonance and the spatial diffusion coefficients. 



Through experimental and theoretical work we show that 
this relation is actually not satisfied by independent mea- 
surements/calculations of the width of the resonance and 
the diffusion coefficients. Instead, we show that the nar- 
row Rayleigh resonance originates from the atomic veloc- 
ity damping, i.e. the width of the resonance is a measure 
of the cooling rate, while spatial diffusion contributes to 
the scattering spectrum with a much broader resonance. 

Consider first the general relation between the width 
of the Rayleigh line and the relaxation rate of the ma- 
terial observable excited in the optical process. In the 
basic setup of pump-probe spectroscopy, an atomic sam- 
ple interacts with two laser fields: a strong pump beam, 
with frequency lo, and a weak probe beam with frequency 
uj + 5. The superimposition of the pump and probe fields 
results in an interference pattern moving with phase ve- 
locity v — S/\Ak\, with Ak the difference between pump 
and probe wavevectors. The atomic sample tends to fol- 
low the interference pattern and a grating of an atomic 
observable (typically density, magnetization or temper- 
ature) is created. However due to the finite response 
time of the atomic medium the material grating is phase- 
shifted with respect to the light interference pattern. 
Therefore the pump beam can be diffracted on the mate- 
rial grating in the direction of the probe, modifying the 
probe transmission. It is then clear that it should be 
possible to derive information about the atomic response 
time from the transmission spectrum. More precisely if 
we assume that only one atomic observable is excited in 
the optical process, and that the time evolution of this 
observable is characterized by a single relaxation rate 7, 
the probe gain spectrum g(S) has then a dispersive line 
shape 
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with peak-to-peak distance 27, as derived in 0- 

Consider now the specific configuration with linearly- 
polarized pump and probe beams, the two polarization 
vectors being parallel. The resulting intensity interfer- 
ence pattern gives rise, via the dipole force, to a grating 
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of the atomic density n of the form 

n = no + ni(exp[— i{8 ■ t — Ak ■ r)] + c.c.) . (2) 

We assume now, following previous work 5], that the 
Rayleigh resonance originates from the scattering on this 
atomic density grating. It follows that the width of the 
Rayleigh line is related to the spatial diffusion coeffi- 
cients. Indeed the relaxation mechanism of a grating of 
atomic density is spatial diffusion: atoms have to move 
to destroy the density grating. More quantitatively, if 
we assume that the dynamics of the atoms in the optical 
lattice is well described by the Fick law 
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where Di (i — x,y,z) is the spatial diffusion coefficient 
in the i direction, we find substituting the expression 
for n in © that the relaxation rate of the atomic 
density, defined by 



is given by 
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Under the assumption that the Rayleigh resonance orig- 
inates from the scattering on the atomic density grating, 
the half-distance peak-to-peak of the Rayleigh line -fR 
is simply equal to the relaxation rate 7_d, and therefore 
measurement of "fu allows the determination of the dif- 
fusion coefficients, as in Refs. H0. The validity of this 
approach will be tested by comparing results for the re- 
laxation rate jrj with measurements of the width of the 
Rayleigh resonance, as presented below. 

In our experiment rubidium atoms are cooled and 
trapped in a three-dimensional (3D) lin_Llin near reso- 
nant optical lattice Q- The periodic structure is deter- 
mined by the interference of four linearly polarized laser 
beams, arranged as in Fig. [2 The angle 29 between 
copropagating lattice beams is equal to 60°. This four- 
beam configuration is the same, except for the value of 
the angle 9, as the one considered in the works of Guibal 
et al p| and Jurczak et al |(|. 

To determine in a direct way the spatial diffusion co- 
efficients of the atoms in the optical lattice, we observe 
the atomic cloud expansion by using a Charge Coupled 
Device (CCD) camera d, U . Since the x and y direc- 
tions are equivalent in our lattice (see Fig. we chose 
to take images in the £Oz plane, where £ is the axis in 
the xOy plane forming an angle of 45° with the x and y 
axis. Correspondingly, we determined the diffusion coef- 
ficients and D z in the £ and z directions, with results 
as in Fig. |21 

These values for the diffusion coefficients are not con- 
sistent with the value of about 10 H/M determined for the 
same configuration by Jurczak et al 6] by polarization- 
selective intensity correlations. As we will show in the 




FIG. 1: Sketch of the experimental setup. 
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FIG. 2: Experimental results for the spatial diffusion coeffi- 
cients in the £ and z directions as functions of the intensity per 
lattice beam II and for different values of the lattice detuning 
A. 



following, this unconsistency derives from the unreliabil- 
ity of the determination of the diffusion coefficients by 
light scattering measurements, as this derivation of the 
diffusion coefficients is based on the assumption that the 
narrow Rayleigh resonance originates from the diffraction 
on an atomic density grating. 

We turn now to the measurements of the width of the 
Rayleigh resonance. The y polarized probe beam is de- 
rived from the lattice beams, with the relative detuning 
controlled with acousto-optical modulators. This probe 
beam is sent along the z axis through the cold atomic 
sample, (Fig. ^) with its frequency scanned around the 
lattice beams' frequency. The probe can interfere with 
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the different lattice beams, which play the role of the 
pump. 

A typical probe transmission spectrum is shown in Fig. 
The lateral resonances have been characterized in 
great detail in past investigations > an d we focus here 
on the resonance at the center of the spectrum (inset of 
Fig. |2J) . To determine whether this Rayleigh resonance 
can be associated with the relaxation mechanism of spa- 
tial diffusion, we made a systematic study of the width of 
the resonance as a function of the interaction parameters 
(lattice- field intensity and detuning). 
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FIG. 3: Probe transmission as a function of the detuning 
between probe and lattice beams, T and To being the intensity 
of the transmitted probe beam with and without the atomic 
cloud. The inset shows a slow scan of the region around zero 
detuning, together with the fit with the function (JSJ (solid 
line). 

The peak-to-peak distance 27^ of the Rayleigh reso- 
nance has been determined by fitting the central part of 
the probe transmission spectrum with the function 




f(5) = ai + a 2 S + 
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Here the linear term in detuning describes the wings of 
the sideband resonances |ll| . The Lorentzian resonance 
arises from the radiation pressure, and has the same 
width of the dispersive line, as discussed in Ref. 
Experimental results for are reported in Fig. 0] as a 
function of the lattice beam intensity, for different values 
of the lattice detuning. 
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FIG. 4: Left: Experimental results for the half of the peak-to- 
peak distance of the Rayleigh resonance. Right: Relaxation 
rate 70 of the atomic density, as calculated from the exper- 
imental data for and D z using Eq. (JHJ. Both quantities 
are plotted as functions of the intensity per lattice beam, for 
different values of the lattice detuning. 



We now describe the determination of 723. In the ex- 
amined configuration the probe beam can interfere si- 
multaneously with all lattice beams. Therefore the situ- 
ation is slightly more complicated than the one analyzed 
previously leading to Eq. [SJ and to derive the link be- 
tween the width of the Rayleigh resonance and the dif- 
fusion coefficients we have to calculate the interference 
pattern between the probe and the lattice beams. By 
using the expression for the lattice-beams electric fields 
for a 3D lin_Llin optical lattice 0] , we easily find that the 
intensity-modulation produced by the probe beam is: 

S\E\ 2 ~ E*E cos(Kx)ex.p{i[(K + -k)z+54]}+c.c. , (7) 

with Eq (E p ) the amplitude of the lattice (probe) field, 
K = k sin 9 and K + = k cos 9. Substituting in Fick's law, 
Eq. Q, the resulting modulation for the atomic density 
we find that the relaxation rate 70, defined via Eq. 
is in the present case 



7_d = D x {ksm9) 2 + D z k 2 {I - cos 9f 



(8) 



This equation is consistent with the relation derived in 
Ref. |5j in the limit of small 9. To determine whether 
the rate of relaxation of the atomic density is equal to 
the width of the Rayleigh resonance, we calculate from 
the values Dj, D z of Fig. [3 the relaxation rate 70 of the 
atomic density, using Eq. ©, with results as in Fig. ^ 
For the same range of interaction parameters the relax- 
ation rate 7d is four orders of magnitude larger than the 
half-distance peak-to-peak 7^ of the Rayleigh resonance. 
We therefore conclude that the Rayleigh resonance does 
not originate from the diffraction on an atomic density 
grating, and therefore measurements of the width of the 
Rayleigh line do not allow the determination of the spa- 
tial diffusion coefficients. 

Our conclusions, based on the presented experimental 
findings, are supported by numerical calculations. We 
consider a J g = 1/2 — > J e = 3/2 atomic transition, 
as customary in numerical analysis of Sisyphus cooling. 
Taking advantage of the symmetry between the x and 
y directions (see Fig. [TJ, we restricted the atomic dy- 
namics to the xOz plane. Through semiclassical Monte 
Carlo calculations [EJ, [L3J , we simulate the dynamics of 
the atoms in the optical lattice. From the atomic trajec- 
tories we determine then the probe transmission spectra 
and the evolution of the atomic mean square displace- 
ments. We calculate the width of the Rayleigh line and 
the spatial diffusion coefficients. From these diffusion co- 
efficients we then derive through Eq. © the relaxation 
rate of the atomic density. The comparison between the 
numerically calculated 7^ and 70, as shown in Fig. [5] 
confirms that the width of the Rayleigh line does not 
correspond to the rate of relaxation of the atomic den- 
sity. 

The final step of our analysis consists in identifying the 
damping process which leads to the phase shift producing 
the Rayleigh scattering. Inspired by previous studies of 
stimulated Rayleigh scattering in corkscrew optical mo- 
lasses |l Ij , we numerically examined the damping process 
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FIG. 5: Numerically calculated relaxation rate of the atomic 
density 70, half-distance peak-to-peak 7k of the Rayleigh line 
and relaxation rates Tt x , Ft z of the atomic temperature in 
the x and z directions. All these quantities are reported as 
functions of the optical pumping rate V' , for different values 
of the light shift per beam A' . Here ui r is the atomic recoil 
frequency. 



of the atomic velocity in the optical lattice and calculated 
the relaxation rates Tj^ , IV, of the atomic temperature 
in the x and z directions, with results as in Fig. |3J It 
appears that the damping rates of the atomic tempera- 
ture not only are of the same order of magnitude of the 
width of the Rayleigh line, but that they also display 
the same linear dependence on the optical pumping rate 
Tq, at fixed light shift per beam A' , i.e. at fixed depth 
of the potential wells. More precisely, neglecting Tx x as 
IV << IV, we find from the data of Fig. [SJthat 



j R = 0.13(±0.04)wr + 0.25(±0.02)r T ^ 



(9) 



This shows that for an optical lattice the width of the 
Rayleigh line is a measure of the cooling rate, a behaviour 
already encountered in corkscrew optical molasses. It is 
then legitimate to investigate the eventual contribution 
of the light scattering on the density grating to the probe 
transmission spectrum. By fitting the broad wings of the 
numerically calculated spectra with a dispersive function, 
we found that the corresponding width is approximately 
equal to the relaxation rate 70 . This shows that the 
scattering on the density grating contributes to the probe 
transmission spectrum with a resonance much broader 
than the narrow line observed at the center of the spectra. 
In other words, the information on the spatial diffusion 
coefficients is contained in the broad wings (ui > 10 MHz) 
of the scattering spectrum, and not in the central narrow 
resonance. 

In summary, in this work we investigated the connec- 
tion between Rayleigh scattering and the atomic dynam- 
ics in dissipative optical lattices. In particular, following 
recent proposals [a, |(| , we studied whether the Rayleigh 
resonance originates from the diffraction on a density 
grating, and is therefore a probe of transport of atoms 
in optical lattices. It turns out that this is not the case: 
the Rayleigh line is instead a measure of the cooling rate, 
while spatial diffusion contributes to the scattering spec- 
trum with a much broader resonance. 
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